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ABSTRACT 

Aims. Massive B-type stars with strong magnetic fields and fast rotation are very rare and provide a mystery for theories of both 
star formation and magnetic field evolution. Only two such stars, called cr Ori E analogs, were previously known. Recently, a team 
involved in APOGEE, one of the Sloan Digital Sky Survey III programs, announced the discovery of two additional rigidly rotating 
magnetosphere stars, HD 23478 and HD 345439. The presence of magnetic fields in these newly discovered cr Ori E analogs was not 
investigated in the past. 

Methods. In the framework of our ESO Large Programme, and one normal ESO programme, we carried out low-resolution PORS 2 
spectropolarimetric observations of HD 23478 and HD 345439. 

Results. From the measurements using hydrogen lines, we discover a rather strong longitudinal magnetic field of the order of up to 
1.5 kG in HD 23478, and up to 1.3 kG using the entire spectrum. The analysis of HD 345439 using four subsequent spectropolarimetric 
subexposures does not reveal the presence of a magnetic field at a significance level of 3cr. On the other hand, the inspection of 
individual subexposures indicates that HD 345439 may host a strong magnetic field, rapidly varying over 88 minutes. A hint at the 
fast rotation of HD 345439 is also given by the behaviour of several metallic and He i lines in the low-resolution FORS 2 spectra, 
showing profile variations already on such a short time scale. 

Key words, stars: early-type — stars: fundamental parameters - stars: individual: HD 23478 and HD 345439 - stars: magnetic field 
- stars: variables: general 


1. Introduction 

The increasing rate of detections of magnetic fields in massive 
stars during the last years generates a strong motivation to study 
correlations between evolutionary state, rotation velocity, and 
surface composition (e.g. Brott et al. I201 II Potter et al. I2012II 
aiming at understanding the origin and the role of magnetic fields 
in massive stars. Recent observations indicate that the presence 
of magnetic fields may be responsible for a wide range of phe¬ 
nomena observed in massive stars, such as chemical peculiarity, 
periodic UV wind-line variability, cyclic variability in Ha and 
He II 44686, excess emission in UV-wind lines, and unusual X- 
ray emission (e.g. Schnerr et al. I20081 Oskinova et al. I201 lak 
The variability is likely due to material that is guided by the 
large-scale magnetic field that corotates with the star (e.g. Groote 
& Hunger [1982l l. To properly understand and model all these 
phenomena, the knowledge of the magnetic field behaviour, i.e. 


* Based on observations obtained in the framework of the ESO Prgs. 
19LD-0255(E) and 094.D-0355(B). 


its configuration, strength, and variability, in massive stars is of 
utmost importance. 

The presence of a rigidly rotating magnetosphere in the 
early B-type stars HD 23478 and HD 345439 was recently dis¬ 
covered in the Apache Point Observatory Galactic Evolution 
Experiment (APOGEE; Eikenberry et al. I20141 I using high- 
resolution (R ~ 22 500) near-infrared H-band spectra. The au¬ 
thors detected in the APOGEE bandpass prominent Brackett 
series emission lines with broad double-horned profiles, which 
are characteristic for the presence of a rigidly rotating magne¬ 
tosphere (RRM), like found in the prototype star cr Ori E (see 
their Eig. 4). The strong magnetic field and fast rotation lead to 
the formation of a centrifugally supported magnetosphere with 
rotationally modulated hydrogen line emission explained within 
the RRM model (e.g. Townsend et al. I2005t . Among cr Ori E, 
HD 23478, and HD 345439, the latter appears to display the 
strongest RRM feature. Eikenberry et al. (I2014I) report that the 
optical spectra of HD 23478 and HD 345439 obtained with the 
Hobby-Eberly-Telescope at a medium resolving power of 18 000 
indicate a classification B3IV for the former star and B1V-B2V 
for the latter. While they determine a v sin i value for HD 23478 
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of about 125 + 20ktns HD345439, in contrast, appears to 
be one of the fastest known rotators among magnetic main- 
sequence stars with v sin i ~ 270 ± 20kms“'. Furthermore, the 
authors report that HD 23478 exhibits a broad Her RRM signa¬ 
ture, while for HD 345439 the Ha profile shows a slight bump 
on the red side. 

In this work, we present spectropolarimetric follow-up ob¬ 
servations of HD 23478 and HD 345439 obtained in the frame¬ 
work of our BOB (B helds in OB stars) collaboration (Morel 
et al. 120141 1 using FORS 2 (FOcal Reducer low dispersion 
Spectrograph) mounted on the 8-m Antu telescope of the VLT. 
We intend to search for the presence of strong magnetic helds 
in both stars, which were implicitly detected through the RRM 
interpretation of their NIR spectra by Eikenberry et al. (120141 1. 


2. Magnetic field measurements using FORS 2 
spectropolarimetry 

Low-resolution FORS 2 observations for HD 23478 and 
HD 345439 were obtained in 2014 on November 17 and June 2, 
respectively. FORS 2 is a multi-mode instrument equipped with 
polarisation analysing optics comprising super-achromatic half¬ 
wave and quarter-wave phase retarder plates, and a Wollaston 
prism with a beam divergence of 22" in standard resolution 
mode. We employed the GRISM 600B and the narrowest 
available slit width of 0'.'4 to obtain a spectral resolving power 
of R ~ 2000. The use of the mosaic detector with a pixel size of 
15 yum allowed us to cover a spectral range from 3250 to 6215 A, 
which includes all Balmer lines apart from Ha, and numerous 
He I lines. For each star, four consecutive observations were 
obtained with pairs of position angles separated by 90°, namely 
[-45 °,-h45°,-h45°,-45°,-45°,-h45°,-h 45°,-45°]. The expo¬ 
sure time for each position angle was 20 sec for HD 23478, and 
600 sec for HD 345439. A peak signal-to-noise ratio (SNR) of 
1900 in the hnal Stokes I spectrum was achieved for HD 23478, 
while for HD 345439 we obtained a peak SNR of ~900. 

The VII spectrum is calculated using: 


V if/ /°-/n i r-r ] \ 


( 1 ) 


where -1-45° and -45° indicate the position angle of the retarder 
waveplate and f° and are the ordinary and extraordinary 
beams, respectively. Null prohles, N, are calculated as pairwise 
differences from all available V prohles. From these, 3cr-outhers 
are identihed and used to clip the V prohles. This removes spu¬ 
rious signals, which mostly come from cosmic rays, and also 
reduces the noise. The mean longitudinal magnetic held, {B^}, is 
measured on the rectihed and clipped spectra based on the rela¬ 
tion 


F _ geff e ^ 
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where V is the Stokes parameter that measures the circular po¬ 
larization, / is the intensity in the unpolarized spectrum, geff is 
the effective Lande factor, e is the electron charge, A is the wave¬ 
length, mg is the electron mass, c is the speed of light, dZ/d/l 
is the wavelength derivative of Stokes /, and (B^) is the mean 
longitudinal (line-of-sight) magnetic held. 

The longitudinal magnetic held is usually measured in two 
ways: using only the hydrogen Balmer lines or using the entire 
spectrum including all available lines. To identify any system¬ 
atic differences that could exist in treating the FORS 2 data by 
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Fig. 1. Graphical output of the iraf/idl “Bonn pipeline” obtained 
from the analysis of the FORS 2 data of HD 23478 consider¬ 
ing the hydrogen lines. The top panel presents the derivative of 
Stokes I. The regions used for the calculation of the magnetic 
held are marked by a thick blue line centered at zero. In the 
middle panel, the top spectrum presents the Stokes I spectrum 
arbitrarily normalised to the highest value, the middle red spec¬ 
trum corresponds to the Stokes V spectrum (in %), while the 
bottom blue spectrum is the N spectrum (in %). The green aster¬ 
isks drawn on the Stokes V and N spectra mark the points which 
have been removed by the sigma clipping. The pale blue strip 
drawn behind the N spectrum shows the uncertainty associated 
with each spectral point. As a double-check of the statistical con¬ 
sistency of the data and uncertainties, the thick green bar on the 
left side of the N spectrum shows the standard deviation of the 
N prohle. In the bottom-left panel, we present the linear ht used 
for the determination of the magnetic held value using Stokes V 
(i.e. (Bz))- The red solid line shows the best ht. From the linear 
ht we obtain (Bz) = -1477 + 95 G. In the bottom-right panel, 
we present the same as in the bottom-left panel, but for the null 
profile (i.e. (Nz)). From the linear ht, we obtain (Nz) - 30+64 G. 


different research groups, the mean longitudinal magnetic held, 
(Bz), was derived in both stars using independent software pack¬ 
ages (one developed in Bonn and the other one in Potsdam). For 
the hrst reduction, we used a suite of IRAF (Tody ll993iPl and 
IDL routines that follow the technique, recipes, and recommen¬ 
dations by Bagnulo et al. (120021120121 l2013lFI . The determina¬ 
tion of the mean longitudinal magnetic held using low-resolution 
FORS spectropolarimetry with the second software package de- 


’ Image Reduction and Analysis Facility (IRAF - 
http://iraf.noao.edu/) is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with the 
National Science Foundation. 

^ More details about the applied data reduction and analysis proce¬ 
dure of FORS spectropolarimetric data will be given in a separate work 
(Fossati et al., in preparation). 
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Table 1. Summary of the results obtained for HD 23478 and HD 345439 from the analysis of the hydrogen lines and of the entire 
spectrum. 


Object 

Subexposures 

Group 

(Bz) [G] 

Hydrogen 
a-B (Vz) [G] 

CTjv 

(Bz> [G] 

All 

ctb 

(Vz> [G] 

crN 

HD 23478 

all 

Bonn 

-1477± 95 

15.5 

30+ 64 

0.5 

-1302+ 59 

21.9 

74+ 49 

1.5 



Potsdam 

-1347+114 

11.8 

50+ 93 

0.5 

-1139+ 84 

13.5 

45+ 70 

0.6 

HD 345439 

all 

Bonn 

4±121 

0.0 

23+111 

0.2 

-2+ 78 

0.0 

36+ 74 

0.5 



Potsdam 

57+153 

0.4 

58+136 

0.4 

-22+116 

0.2 

64+103 

0.6 


1+2 

Bonn 

649+169 

3.8 

-117+157 

0.7 

389+105 

3.7 

104+101 

1.0 



Potsdam 

638+183 

3.5 

-197+188 

1.0 

507+134 

3.8 

-93+143 

0.6 


3+4 

Bonn 

-618+171 

3.6 

188+154 

1.2 

-375+109 

3.4 

18+103 

0.2 



Potsdam 

-574+193 

3.0 

281+192 

1.5 

-501+146 

3.4 

207+145 

1.4 


1 

Bonn 

435+244 

1.8 



418+159 

2.6 





Potsdam 

414+282 

1.5 



436+212 

2.1 




2 

Bonn 

631+253 

2.5 



252+148 

1.7 





Potsdam 

789+246 

3.2 



565+188 

3.0 
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Bonn 

-509+242 

2.1 



-402+160 

2.5 





Potsdam 

-303+282 

1.1 



-298+212 

1.4 
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Bonn 

-1181+328 

3.6 



-534+166 

3.2 





Potsdam 

-840+262 

3.2 



-689+198 

3.5 




veloped in Potsdam is described by Hubrig et al. (120141120151 
and by Scholler et al. (in preparation). In general, the measure¬ 
ment method is the same for both pipelines. A few minor dif¬ 
ferences refer to the clipping procedure, rectification, and the 
choice of the wavelength regions including the hydrogen lines. 
Furthermore, Monte Carlo bootstrapping tests are carried out in 
the second software package (e.g. Rivinius et al. 120101) . In these 
tests, we generate 250000 statistical variations of the original 
dataset by the bootstrapping technique, and analyze the resulting 
distribution P((Bz)) of the regression results. Mean and standard 
deviation of this distribution are identified with the most likely 
mean longitudinal magnetic field and its 1 cr error, respectively. 
The main advantage of this method is that it provides an inde¬ 
pendent error estimate. The measurement uncertainties for both 
stars obtained before and after Monte Carlo bootstrapping tests 
were found to be in close agreement, indicating the robustness 
of the measurement method. 

The results of the magnetic field measurements carried out 
using independent reduction and measurement packages are 
listed in Table [T] For the first software package developed 
in Bonn, the uncertainty on the {B^} and {N^} values is the 
;if^-scaled uncertainty (Bagnulo et al. 120121 ). while the uncer¬ 
tainty determined by the second software package developed in 
Potsdam involves in addition the bootstrapping results. In each 
column, the numbers ctb and cr^ after the (B^) and (Nz) values 
refer to the sigma detection. 

A strong magnetic field at a significance level of more than 
lOcr was achieved for HD 23478 using the hydrogen lines and 
the entire spectrum. Using the first software package we obtain 
{Bz} - -1302 + 59G for the whole spectrum and {Bz)hyd = 
-1477 + 95 G, while using the second software package we mea¬ 
sure (Bz) = -1139 + 84 G and (Bz>hyd = -1347 + 114 G. For 
both reduction packages the magnetic field values obtained for 
each pair of position angles separated by 90° show negative po¬ 
larity, are all consistent with each other, and are of similar or¬ 
der. No significant fields were detected in the null spectra cal¬ 
culated by combining the subexposures in such a way that the 
intrinsic source polarization cancels out, yielding a diagnostic 
null N spectrum. Figure [Upresents the main graphical output of 
the “Bonn pipeline” showing the clear magnetic field detection 
obtained form the analysis of the hydrogen lines of HD 23478. 



Fig. 2. Variation of the (Bz) value obtained for HD 345439 from 
the analysis of the hydrogen lines (black asterisks) and of the 
whole spectrum (blue rhombs). Between each point there is a 
time span of about 22 min. 

In the analysis of HD 345439, using four consecutive obser¬ 
vations obtained with pairs of position angles separated by 90°, 
we obtained a clear non-detection, which is presented in the sec¬ 
ond line in Table[T] However, after splitting the observations into 
two data sets, i.e. using the first two pairs and the second two 
pairs consisting of observations at the retarder waveplate posi¬ 
tions [-45°, -)-45°, -1-45°, -45°], we obtain 3.0 to 3.8cr detections 
shown in the third and fourth lines of the table. The fact that 
the first and second data-sets led to opposite (Bz) values, though 
almost equal in modulus, led to a non-detection in the measure¬ 
ments using the whole available data set. This behaviour might 
be suggestive of a very fast rotating star with a magnetic field 
axis inclined relative to the rotation axis. Given the apparent 
fast and large variation of the value of the longitudinal mag¬ 
netic field, we then decided to measure the magnetic field using 
each individual pair obtained at the retarder waveplate positions 
[-45°, -)-45°] and separated by roughly 22 minutes. The results 
are presented in the last four lines of Table [1] In this case, the 
null profile can not be extracted, since the measurements were 
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Fig. 3. Left-side panel: Stokes V spectrum of HD 345439 in the vicinity of Hy and Hei /14471 obtained from the combination 
of the first two pairs of subexposures. The shape of the Zeeman features indicates the presence of a longitudinal magnetic field 
of positive polarity. Right-side panel: Stokes V spectrum of HD 345439 in the vicinity of Hj3 and He i A4922 lines obtained from 
the combination of the last two pairs of subexposures. The shape of the Zeeman features indicates the presence of a longitudinal 
magnetic field of negative polarity. 


done using observations at just two position angles. Note also 
that no sigma clipping was adopted for this specific analysis, be¬ 
cause the sigma clipping is applied on the basis of the null pro¬ 
file. Figure|2]shows the variation of the (B^) value obtained from 
the analysis of the hydrogen lines and of the entire spectrum us¬ 
ing the first software package developed in Bonn. We discover 
that the mean longitudinal magnetic field measured on hydrogen 
lines is changing from about -1-500 G measured in the first pair 
of subexposures to about -1200G measured in the last pair of 
subexposures. Using the second software package developed in 
Potsdam we obtained similar magnetic field behaviour, but with 
somewhat different absolute values of the field. The detections 
are achieved at a significance level of 3 to 3.6cr in two subex¬ 
posures. In Figure 12 we present Stokes V spectra obtained for 
the combinations of the first two pairs of subexposures and for 
the last two pairs of subexposures in the vicinity of the hydrogen 
and He I lines. The different shape of the Zeeman features indi¬ 
cates that the magnetic field polarity has changed after the the 
second subexposure. We note that all tests we performed on the 
data did not give an indication of the presence of relevant sys¬ 
tematic noise. Nevertheless, it is important to obtain in the future 
more spectropolarimetric data to confirm this behaviour. 

An indication for the fast rotation of HD 345439 is also given 
by the behaviour of several metallic and He I lines in the low- 
resolution FORS 2 spectra, showing profile variations already 
on a time scale of 88 min. In Fig. |4] we present the behaviour 
of the Si III 44553 and the He i 44388 and He 1 44471 lines in the 
FORS 2 Stokes I spectra obtained for each pair of subexposures 
recorded with a time lapse of 22min. The profile of the Sim 
line becomes asymmetric in the last pair of subexposures, i.e. at 
the time when we observe negative field polarity. We note, how¬ 
ever, that the level of variability of the Si iii 44553 line profile is 
not significantly higher than that in neighbouring continuum re¬ 
gions. On the other hand, similar asymmetric profiles of Si lines 
are detected in our FEROS (Fiber-fed Extended Range Optical 
Spectrograph) spectrum of cr Ori E obtained in the framework of 
the ESO programme 091 .C-0713(A). As for the behaviour of the 


He I 44388 and Hei 44471 lines, their intensity is lower at the 
phase when the polarity of the field is positive and it becomes 
higher when the field is negative. Such a fast spectral and lon¬ 
gitudinal magnetic field variability suggests that the rotation pe¬ 
riod should be very short, probably less than one day. The rapid 
change in the strength of the mean longitudinal magnetic field 
by about 1 kG over 88 min in HD 345439 is reminiscent of the 
mean longitudinal magnetic field behaviour in the fast rotating 
cr Ori E analog HR 7355 with a rotation period of 0.52 d. This 
star shows a field variation of about 1 kG within a phase inter¬ 
val of 0.1, corresponding to 75 min (see Eig. 3 in Rivinius et al. 
120131) . Therefore, it is not impossible that the rotation period of 
HD 345439 is also significantly shorter than one day. 

No spectral variability was detected in the eight low- 
resolution EORS2 spectra for HD 23478 taken within 
15 min. This star was intensively studied photometrically 
by Jerzykiewicz (119931 ). who suggested a photometric/rotation 
period of 1.05 d and, using a radius estimated from the spectral 
type, obtained 240 km s ' for the equatorial rotation velocity. 


3. Discussion 

Our spectropolarimetric follow-up of the early B-type stars 
HD 23478 and HD 345439, with recently discovered rigidly ro¬ 
tating magnetospheres, using EORS 2 observations, reveals the 
presence of a strong longitudinal magnetic field of the order 
of 1.5 kG in HD 23478. We also found an indication for the 
presence of a rapidly changing magnetic field in the fast ro¬ 
tator HD 345439. Obviously, additional follow-up spectropo¬ 
larimetric monitoring of the behaviour of the magnetic field 
in HD 23478 and HD 345439 is urgently needed to confirm 
and characterize the magnetic field topology of both stars. 
Knowledge of the structure of the magnetic field is of crucial 
importance for the understanding of the nature of these stars be¬ 
longing to the group of cr Ori E analogs. 

This type of stars is extremely rare: apart from crOriE (e.g. 
Groote & Hunger 119971) . only two other stars, HR 5907 and 
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Fig. 4. Line profile variability of Si and He in the FORS 2 Stokes 
I spectra of HD 345439 obtained for each pair of subexposures, 
i.e. with a time lapse of 22 min. Top panel presents the behaviour 
of the Si III 44553 line. The mean profile is indicated by the red 
dashed line. The numbers close to the plotted line prohles re¬ 
fer to the time of the observation calculated as MJD-56810. The 
vertical bar on the left side shows the noise level. In the bottom 
panel we present overplotted profiles of the He i 44388 and He i 
44471 lines. Their intensity is lower at the phase when the polar¬ 
ity of the field is positive and it becomes higher when the field 
is negative. The error bar is too small to be presented in these 
plots. 


HR 7355, are known to belong to this group called the crOriE 
analogs (Grunhut et al. 120121 Rivinius et al. 120131 1. Their very 
strong magnetic fields and extraordinary fast rotation provide a 
mystery for theories of both star formation and magnetic held 
evolution. Both stars have similar dipole strengths of the order of 
10 kG and similar very short half-day rotation periods. However, 
the spindown timescale via magnetic braking is expected to be 
much shorter than their estimated ages (e.g., Rivinius et al. 120131 
Mikulasek et al. 120101) . 

The emission variability in hydrogen lines in the cr OriE 
analogs suggests that gas is trapped in magnetospheric clouds. 
The high strength of the magnetic helds controlling the stellar 
winds and fast rotation makes these type of stars a test case for 
models of X-ray production. Only low resolution X-ray spectra 
were so far measured for three members of this group, cr Ori E, 
HR5907, and HR7355 (e.g.. Petit et al. l20T3l Naze et al. l20T4l) . 
We retrieved and analyzed these archival XMM-Newton and 
Chandra observations using most recent calibrations. Our anal¬ 
ysis of the X-ray spectra of HR 5907 and HR 7355 shows that 
these stars are even harder X-ray sources than cr OriE; e.g., the 
emission measure weighted temperature of the hot plasma in 
HR 5907 is 80 MK, and in HR 7355 it is 60 MK. Given that the 
winds in magnetic B-type dwarf stars are weak, with velocities 
not exceeding a few hundred kms * (Oskinova et al. 1201 IbT l. it 
is very difficult to explain such high temperatures in the frame¬ 
work of the magnetically confined wind shock model (Babel & 


Montmerle 1 1 9971) . On the other hand, non-thermal emission can 
be expected from stars with strong magnetic helds (Leto et al. 
120061) . Our analysis of archival data showed that the presence of 
non-thermal radiation in these objects is very likely (Oskinova 
et al., in preparation). New, already approved XMM-Newton ob¬ 
servations of HD 345439 will probe whether hard non-thermal 
X-rays are present in this object and shed new light on the mech¬ 
anisms of X-ray generation in strongly magnetic fast rotators. 
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